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Abstract. We develop the white dwarf luminosity func- 
tion (LF) of the old open cluster NGC 188 in order to 
determine a lower limit to the age of the cluster by using 
the faint end of the cooling sequence. 

To produce an extensive sequence of the cooling white 
dwarfs we imaged four contiguous HST-WFPC2 fields in 
the center of the cluster in the F555W and F814W filters. 
After imposing selection criteria on the detected objects 
we found a white dwarf cooling sequence (down to V ~ 
26.5) including 28 candidate white dwarfs in the cluster. 
The exposures are not deep enough to reach the end of this 
sequence, but the results of our analysis allow us to estab- 
lish a lower limit to the age of the cluster independently 
of the isochrone fit to the cluster turnoff. 

The most ancient white dwarfs found are ~ 4 Gyr old, 
an age that is set solely by the photometric limit of our 
data. Classical methods provide an estimate of ~ 7 Gyr 
(Sarajedini et al., 1999). 

Key words: Open clusters and associations: individual 
(NGC 188)- stars: white dwarfs - stars: luminosity func- 
tion 



1. Introduction 

Stellar evolution theory predicts that all single stars hav- 
ing a Main Sequence (MS) mass lower than ~ 8 M Q end 
their lives as white dwarfs (WDs). In a star cluster, the 
white dwarf population carries key information for ad- 
dressing a number of astrophysical questions. In partic- 
ular, since the white dwarfs cool at a predictable rate, 
their LF can be used to constrain the cluster age. This 
age estimate can provide an independent test of the more 
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traditional age determinations based on models and ob- 
servations of the main-sequence turnoff (MSTO). 

The white dwarf cooling sequence in the color- 
magnitude diagram, in combination with models can be 
used to estimate the initial mass - final mass relation- 
ship between the stellar progenitors and the white dwarfs 
over the range of progenitor mass from near 8M Q to the 
current-day turnoff mass of a cluster. This allows a de- 
termination of the amount of mass lost from stars during 
their evolution through stellar winds and planetary nebula 
ejection. 

Ideally, to employ white dwarfs as age estimators, we 
need to observe a large enough sample in any cluster to 
define the WDLF sufficiently well to pinpoint the luminos- 
ity of the turndown with a precision equivalent to a Gyr 
of cooling time. With smaller samples, it is still possible 
to make some progress on the age questions. For example, 
even a single good candidate WD in a cluster at magni- 
tudes fainter than the faintest WDs predicted from the 
MSTO age and cooling theory can bring into question the 
vailidity of the WD or MSTO-based ages. 

Clearly the most interesting case for testing MSTO- 
based ages is for the globular clusters. However, because 
of the combination of large age and large distance for 
even the nearest of the Galactic globulars, it is not triv- 
ial to reach the end of the cooling sequence (e.g. Richer 
et ah, 1995, 1997, 2002; Hansen et al. 2002). Old open 
clusters provide an interesting alternative for WD-cooling 
age studies. There are several populous clusters which are 
both nearer and younger than the nearest globulars and 
for which the expected WDLF turndown is 3 magnitudes 
or more brighter. For open clusters older than ~ 4 Gyr, 
the main-sequence turnoff stars are very similar in struc- 
ture to globular cluster turnoff stars, most importantly 
in that they have radiative cores, and comparison of nu- 
clear and cooling ages for these old open clusters are very 
relevant the issue of globular cluster MSTO age tests. A 
significant difficulty with open clusters is that even the 
most populous are sparse compared to most globular clus- 
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ters. This leads to the problem that the lowest luminos- 
ity WDs are typically dramatically outnumbered by faint 
blue galaxies in the cluster fields. Attempts to statistically 
derive a WDLF by subtracting counts from nearby con- 
trol fields are compromised by small errors in control-field 
counts and cosmic dispersion in faint blue galaxy surface 
density. Very accurate star-galaxy separation is therefore 
crucial for open cluster WD studies. Given the compact 
size of many faint blue galaxies, Hubble Space Telescope 
(HST) imaging gives a tremendous advantage. Ultimately 
another powerful technique for isolating cluster-member 
WDs from the compact blue galaxy and field WD popu- 
lations may be measuring proper motions. 

NGC 188 was considered for a long time to be the old- 
est observable open cluster in the Galaxy. Almost 40 years 
ago Sandage gave an estimate fo r the age of this cluster o f 
10 Gyr (Sandage, 1962; see also Eggcn fc Sandage, 1969| ). 
Recent papers revise the early age estimates downward, 
and more modern values give an age near 7 Gyr (see 
Sarajedini et al., 1999). Other open clusters have now 
been found to be older than NGC 188 (NGC 6791 and 
Berkeley 17 with ages of 8 and 12 Gyr respectively, see 
Phelps et al., 1994J ) however, NGC 188 is still one of the 
most studied since it is both nearer and less obscured than 
some of the other old open clusters. 

Sects. H and [?] present the data, the reduction pro- 
cedures and the color magnitude diagrams. The resulting 
LFs together with final results are presented and discussed 
in Sects. | and |. 

2. Observations and data reduction 

Deep HST WFPC2 images in F555W and F814W were 
obtained around the centre of NGC 188 in the Cycle 7 pro- 
gram GO 7371. The observations consisted of four fields, 
fl, f2, f3 and f4 located as shown in Table [l], in such a way 
that WF4 of the field f4 overlaps WF2 of the field fl and 
WF4 of the field f3 overlaps WF2 of the field f2. Table 
also lists the corresponding exposure times together with 
the number of observed frames for each filter. 

Table 1. Journal of observations 
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Corrections to the raw data for bias, dark and flat- 
fielding were performed using the standard HST pipeline. 
Before reduction, we masked bad pixels and corrected the 
images for vignetting by using a mask and a vignetting 
mask respectively as built for HST data (supplied by P. 
Stetson). We also multiplied our science exposures by a 



pixel-area map to correct for geometrical distortion and to 
restore the integrity of the flux measurements. Photome- 
try was measured using ALLFRAME (Stetson-version 3, 
1997) with a quadratically varying point spread function 
(PSF). 

The images did not have enough bright unsaturated 
stars to create proper PSFs, so PSFs for each chip were 
constructed using WFPC2 images of the globular cluster 
ui Cen images (Stetson P. B. 1997, private communica- 
tion). Since there were virtually no stars in the PC frames, 
these were not included in the reduction and analysis. The 
data reduction process was completed by doing the charge 
transfer efficiency (CTE) correction on the data-set by us- 
ing an algorithm supplied by P. B. Stetson (priv. comm.). 

Finally the instrumental magnitudes were trans- 
formed into the standard Johnson system following 
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To separate stars from other objects at faint magni- 
tudes below V ~ 24.5, we applied SExtractor to the images 
(Bcrtin et al., 1996). This program uses a neural network 
to distinguish between faint point sources from extended 
objects. It assignes a stellarity index c to the objects, rang- 
ing from (galaxies) to 1 (stars). Figure [l] shows the re- 
sults of the SExtractor morphological classification versus 
the V magnitude for all the objects in our sample. Look- 
ing at the figure there is a clear separation between the 
objects with stellarity index around 1 and those with stel- 
larity index around and there is very little in between. 

Careful examination of each object by eye showed that 
virtually all the objects for which < c < 0.4 are classi- 
fiables as galaxies, hot and warm pixels, ghost images or 
part of diffraction spikes; all the objects for which 0.4 < 
c < 0.8 are dubious objects; all the objects for which c > 
0.8 may be considered stars. 

Using these criteria, we established a catalog of 157 
stars in all the chips and all the fields studied selected 
from a total number of 1342 objects. 



3. Data analysis 

3.1. The colour-magnitude diagram and the WD cooling 
sequence 

The V, V — I CMD for all stars (as indicated from the 
stellarity index) measured in the four fields is shown in 
Fig. | 

There is a fairly well-defined main sequence extending 
almost 6 magnitudes from V ~ 19 down to V ~ 25. The 
solid line overlaid on our data is the ridge line obtained by 
using a sample of ground-based data (von Hippel et al., 
1998). 

In the lower left of the CMD there are 28 blue objects 
ranging from V ~ 23 with V-I ~ to V ~ 26.5 with 
V-I ~ 1. An enlargement of this region of the CMD is 
shown in Fig. ||; for each star we also plotted the associated 
photometric error as computed by ALLFRAME. 
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Fig. 1. SExtractor stellarity index vs the magnitude V 
for objects in our sample. The horizontal dashed lines at 
c = 0.4 and c = 0.8 show the thresholds we used to dis- 
tinguish between real stars and other objects (i.e. cosmic 
rays, image defects, galaxies). 



Note that in this figure all objects appear to be stars 
but, for now, we cannot discriminate between those be- 
longing to the cluster and those belonging to the Galactic 
field. We expect a contamination of ~ 8 Galactic white 
dwarfs at the magnitude and colour of interest for our 
analysis (see Mendez & Minniti, 2000). The correction for 
the field contamination will be addressed in the following 
sections. 

We can use the CMD to estimate the expected number 
of WDs in the cluster, starting from the number of the 
MS stars. To carry out this consistency check, we have 
found the slope of the mass function of the cluster by 
using the number of MS stars in the V magnitude range 
20 - 24, obtaining a value of 1.72 (Salpeter is 2.35). We 
then estimated the mass of the stars at the TO under the 
hypothesis that the cluster age is 4 - 5 Gyr, obtaining 
respectively ~ 1.23 Mq and 1.15 M©. By looking at the 
coolest and faintest WDs inside the range in magnitude V 
(23 - 26), and adopting a WD mass of 0.6 M©, we obtained 
only 3 and 4 expected WDs in the cluster (respectively 
for 4 and 5 Gyr), a very small number if compared with 
the observed number in the same range in magnitude (23 
WDs). 

We emphasize that this result is a consequence of 
the assumptions on the poorly determined mass func- 
tion slope. For example, for a slope of 1.05, this number 
changes to 13 or 22 for cluster ages of 4 or 5 Gyr re- 
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Fig. 2. V — I colour-magnitude diagram for the objects 
in the four observed fields. The solid line overlayed on the 
data is the ridge line obtained by using the data sample 
from von Hippcl et al. (1998). 



spectively. If we then are to believe that the most of the 
objects in the observed cooling sequence are in fact WDs, 
we need to assume that our sample is somehow deficient 
in MS stars. This is not unexpected for NGC 188 because, 
as suggested by many authors (Sarajedini et al., 1999, von 
Hippel et al., 1998), dynamical processes and mass segre- 
gation are important in the history of this cluster which 
has the rather short relaxation time of t r h = 6.4 x 10 7 yr 
(McClure et al., 1977). 

3.2. Completeness correction 

Figure || also indicates the completeness in the counts 
computed in bins of 0.5 magnitude in V . Incompleteness 
corrections were determined by carrying out artificial star 
tests on both sets of the V and I frames, so as to be able 
to estimate the overall completeness of our final colour 
magnitude diagram in the region of the candidate white 
dwarfs. To accomplish this we first applied artificial star 
tests to the I images. Artificial stars were added to ran- 
dom positions in the images over a range of 3 magnitudes 
in the region of the CMD of interest for the candidate 
white dwarfs. We created 5 artificial images for each orig- 
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We have estimated the correction to our WD candi- 
date counts, by assuming that all resolved galaxies have 
already been eliminated from our data-sample. We then 
used the number of unresolved blue objects identified by 
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Fig. 3. V — I colour-magnitude diagram for all possible 
white dwarf candidates in our sample. Objects were se- 
lected based on their stellarity index (c > 0.8) and, at 
this stage, there is no correction for field contamination. 
For each object we have also plotted the associated error 
based on the frame-to-frame scatter in the measurements. 
Dotted lines (and the percentages) on the left indicate the 
mean completeness for each bin at the respective V mag- 
nitude. 



inal image each with 100 artificial stars. We emphasize 
that crowding was not a problem in our images. 

We then added an equal number of stars at the same 
position of the artificial / stars to the F-band frames and 
with a magnitude such that they would fall along the white 
dwarf cooling sequence. All the V and / frames with artifi- 
cial stars were then subjected to the same analysis adopted 
for the original frames. 

An artificial star was considered recovered if <5X, 8Y < 
1.5 pixel and <5mag < 0.3 mag. The ratio N lec /N s i m = $, 
taken in 0.5 magnitude bins gives the completeness in that 
bin for the image considered, where N rcc is the number of 
the recovered stars and N s j m is the number of simulated 
stars. 

3.3. Accounting for the field white dwarfs 

In addition to correction for photometric incompleteness 
it is necessary to account for the contamination of the field 
stars. The majority of the field star contamination is due 
to main-sequence stars from the Galactic disk together 
with a few main-sequence stars from the thick disk and 
the halo. The majority of these objects have much redder 
colors than our WD candidates in NGC 188. 



Mendcz fc Minniti (2000)| in the Hubble Deep Field North 
in an area ~ 4 times smaller than that covered by our 
fields. 

Table ^ lists the counts for 0.5 magnitude- wide bins in 
V, the completeness factor and the counts corrected for 
this factor. The last two columns report the number of 
unresolved blue objects for each bin and our counts after 
the correction for completeness and for the contribution 
of the field. 

The absolute V magnitudes of Col. 2 have been obtained 
adopting (m-M) v = 11.44, E B -v=0-09 (Sarajedini et al, 
1999). 
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Table 2. Luminosity function for the candidate white 
dwarfs in our sample in the V band in 0.5 mag bins. N is 
the actual number of stars observed, $ is the completeness 
correction factor in each bin, N c is the number of stars 
after the corrections for incompleteness have been applied, 
Nfidd is the number of unresolved blue objects as found 
in Mendez & Minniti (2000) for the magnitude and the 
colour of our sample, N s is the number of stars after the 
incompleteness correction and the field star subtraction, 
i.e. the cluster white dwarf luminosity function. 



4. Discussion 

To estimate the WD cooling age for NGC 188, our first ap- 
proach is to compare the observed WD CMD with the the- 
oretical WD isochrones corresponding to different cluster 
ages. Figure [| shows the theoretical WD isochrones corre- 
sponding to cluster ages of 1, 2, 3, 4, 5, 6 and 7 Gyr (solid 
lines), superimposed on the observed NGC 188 WD cool- 
ing sequence. In the figure we also plot the theoretical cool- 
ing sequences for 0.5, 0.7 and 0.9 M Q hydrogen-rich WD 
(dotted lines) to indicate the mass range of the WDs in 
the cluster. In particular most of the WDs cover the mass 
range included in our models while a small fraction seem 
to have a mass larger than 0.9 M . Isochrones and the- 
oretical cooling sequences have been developed using the 
WD cooling models of Hansen (1998, 1999) as described 
in Richer et al. (2000) and have been transformed to ap- 
parent magnitudes and colors appropriate for NGC 188 
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adopting E(V-I) = 1.21 E(B-V) (based on equations 3a 
and 3b of Cardelli et al., 1989 applied to the F814W fil- 
ter) and (m-M)v = 11.44. 




Fig. 4. Comparison of our NGC 188 cooling sequence 
with theoretical isochrones (solid lines) and hydrogen-rich 
WD cooling sequences (dashed lines) obtained by using 
models from Hansen (1998, 1999). 

Since the population of WDs accumulate at lower lu- 
minosities as the age of the cluster increases, we can fix a 
limit to the age of the cluster by simply looking at the loca- 
tion of the WDs with respect to the theoretical isochrones. 
Figure || seems to show that for NGC 188 this limit is ~ 
6-7 Gyr. However this method should be used with cau- 
tion. In fact, when we compare the observed WD CMD 
and the theoretical isochrones we do not take into account 
two important corrections that must be applied to the 
data before earring out any quantitative conclusion, i.e., 
the correction for incompleteness and the subtraction of 
the field contamination. These corrections may be negligi- 
ble when the sample is composed of many stars but they 
become very important when we are dealing with such few 
objects. 

The second approach we can follow to derive an esti- 
mate of the age of NGC 188 is that of studying the ob- 
served WDLF. 

As already mentioned above, the population of WDs 
tend to accumulate at lower luminosities as the age of the 
cluster increases. Hence the corresponding white dwarf LF 
will show a peak whose location in magnitude is a function 



of the cluster age. In particular, the peak of the white 
dwarf LF moves towards fainter luminosities as the cluster 
ages (see Richer et al. 2000). 

Figure [| (solid line) shows the differential observed 
WDLF after corrections for incompleteness and field con- 
tamination have been applied (see table ||). In the figure 
error bars reflect the total error associated with each bin 
and include both the counting uncertainty and the uncer- 
tainty due to the correction for the incompleteness: 
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where N is the number of observed WDs; $ is the in- 
completeness factor and N s ; m is the number of simulated 
stars. 

The uncertainty in the number of objects for each bin 
has been estimated assuming that the counting uncer- 
tainies are derived from a Poisson distribution, and that 
the uncertainties in determining $ are derived from a bi- 
nomial distribution as shown in Bolte (1989). 
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Dashed lines in the four panels of the figure refer to 4 
synthetic LFs corresponding to 4, 5, 6 and 7 Gyr respec- 
tively. They have been computed by adopting a Salpeter 
initial-mass function (Salpeter 1954, n(m) oc m~ Q ;a = 
2.35) as described in Richer et al. (2000) and have been 
normalized to the total number of WDs brighter than 
14.75, which is the magnitude of the faintest observed WD 
in the observed LF. 

Figure || shows that the peak of the theoretical LF is 
located at magnitudes fainter than the last observed data 
point (My=14.75) only in the 6 and 7 Gyr cases. Since it 
is unlikely that our data reach the true end of the cooling 
sequence, we suggest this is an indication that this cluster 
is at least as old as 5 Gyr. 

5. Conclusions 

Although we have identified a clear population of objects 
in NGC 188 with the photometric properties expected of 
cluster WDs, it is apparent from Fig. || that we have not 
reached the turnover in the WDLF. The comparison with 
the LF models suggests that we can place a lower limit of 
5 Gyr on the age of the cluster. 

These observations could be improved significantly 
with a modest effort, with HST and the ACS. A fainter 
photometric limit together with an expanded field cover- 
age would allow a large sample of WDs to be identified in 
the NGC 188 field. NGC 188 stars have a proper motion 
of 2.3 milli-arcseconds with respect to the field and with a 
followup program, a subsample of bona-fide cluster WDs 
could be identified on the basis of their proper motions. 



6 



G. Andreuzzi et al.: White Dwarfs in NGC 188 




Richer, H., B., Hansen, B., Limongi, M. et al, 2000, ApJ 529, 
318 



Richer, H., B., B rewer, J., Fahlman, G. G. et al., 2002, |astro- 



ph/0205086 



Salpeter, E., E., 1954, ApJ 121, 161 
Sandage, A., 1962, ApJ 135,333 

Sarajedini A., von Hippel, T., Kozhurina-Platais V., & Demar- 

que, P., 1999 AJ 118, 2894 
von Hippel, T., & Sarajedini, A., 1998 AJ, 116, 1789 
von Hippel, T., & Gilmore, G., 2000 AJ 120, 1384 



11 12 13 14 15 11 12 13 14 15 



Fig. 5. Differential NGC 188 white dwarf luminosity func- 
tion (solid line) , compared with synthetic luminosity func- 
tions of ages 4, 5, 6 and 7 Gyr (dashed lines). 
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